The Electromagnetic Self-Energy Contribution to M v — M n 
and the Isovector Nucleon Magnetic Polarizability 
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We update the determination of the isovector nucleon electromagnetic self-energy, valid to leading 
order in QED. A technical oversight in the literature concerning the elastic contribution to Cotting- 
ham's formula is corrected and modern knowledge of the structure functions is used to precisely 
determine the inelastic contribution. We find 5M~' n = 1.30(03) (47) MeV. The largest uncertainty 
arises from a subtraction term required in the dispersive analysis, which can be related to the isovec- 
tor magnetic polarizability. With plausible model assumptions, we can combine our calculation with 
additional input from lattice QCD to constrain this polarizability as: j3 p - n — —0.87(85) x 10 _4 fm 3 . 

PACS numbers: 13.40.Dk, 13.40.Ks, 13.60.Fz, 14.20.Dh 



Given only electrostatic forces, one would predict that 
the proton is more massive than the neutron but the 
opposite actually occurs [l|-[3 : 



M„ - M p = 1.29333217(42) MeV . 



(1) 



Before we knew of quarks and gluons there were many 
attempts to explain this contradiction, see Ref. |4j for a 
review. We now know there are two sources of isospin 
breaking in the standard model, the masses of the up 
and down quarks as well as the electromagnetic interac- 
tions between quarks governed by the charge operator. 
The effects of the mass difference between down and up 
quarks are larger and of the opposite sign than those of 
electromagnetic effects, see the reviews 0t3- The net 
result of the quark mass difference and electromagnetic 
effects is well known, Eq. ([1]), but our ability to disen- 
tangle the contributions from these two sources remains 
poorly constrained. 

In contrast, lattice QCD calculations have matured sig- 
nificantly There are now calculations performed with the 
light quark masses at or near their physical values [9-[l2| , 
reproducing the ground state hadron spectrum within 
a few percent. These advances have allowed for calcu- 
lations to begin including explicit isospin breaking ef- 
fects from both the quark masses [13| - K17| and electro- 
magnetism [TJI, [l§ - l2ll ]. While the lattice calculations 
of rrid — m u effects are robust, the contributions from 
electromagnctism are less mature and suffer from larger 
systcmatics, due in large part to the disparity between 
the photon mass and a typical hadronic scale. Improved 
knowledge of rrid — in u and its effects in nucleons will en- 
hance the ability to use effective field theory to compute 
a variety of isospin- violating (charge asymmetric) effects 
in nuclear reactions [7|, I22l427j | . 

An application [28[ of the Cottingham sum rule |29| > 
which relates the electromagnetic self-energy of the nu- 
cleon to measured elastic and inelastic cross sections, 
gives the result SM^_ n = 0.76 ± 0.30 MeV. Given the 



high present interest in the precise value of <5M7 n 



and 



its many possible implications, it is worthwhile to revisit 
this result. Many high quality electron scattering exper- 
iments have been performed since 1975 and there have 
also been theoretical advances. The central aim of this 
work is to provide a modern, robust evaluation of SMp_ n . 
We will show the precision of this effort is severely lim- 
ited by our knowledge of the required subtraction func- 
tion. Given plausible model assumptions, this limitation 
is translated into our knowledge of the isovector nucleon 
magnetic polarizability, /3 p _ n = 0%, — ft^, for which even 
the sign is presently unknown [301 ]. 

Cottingham 's sum rule- In perturbation theory, the 
electromagnetic self-energy of the nucleon, <5M 7 , can be 
related to the spin averaged forward Compton scattering 
tensor 

T ^ = \llj d ^ e*«(HT WOJ,(0)} b*> . (2) 
integrated with the photon propagator over space-time 
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2AI (2tt) 3 



A- 



(3) 



where we work in the nucleon rest frame p^ — (M, 0), 
a = e 2 /4-7r and the subscript R implies the integral has 
been renormalized. Performing a Wick rotation of the 
integration contour to imaginary photon energy, the nu- 
cleon self-energy can be related to the structure functions 
arising from the scattering of space-like photons through 
dispersion theory, giving rise to what is known as Cot- 
tingham's formula (the Cottingham sum rule) [29|,[3l|]. In 
principle, this allows the integral in Eq. ([3]) to be com- 
puted in a model independent fashion with input from ex- 
perimental data. There are a few issues which complicate 
the realization of this method: a subtracted dispersive 
analysis is required introducing an unknown subtraction 
function [321, (33[; the integral in Eq. ([3]) diverges loga- 
rithmically in the ultra-violet region and requires renor- 
malization [34|. We review these issues briefly. 
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Lorentz invariancc significantly constrains the form of 
T M „, for which there are two common parameterizations, 



T^(p,q) = -DftTtiv, -q 2 ) + D%T 2 (v, -q 2 ) 



d$q 2 tl(v,-q 2 )-d$q 2 t 2 (v,-q 2 ) 



'(2 



(4a) 
(4b) 



where p ■ q = Mv and 

d (i) - D {i) =a _ Mil 

f 



Ijp. \rprv 2 W" ' f<"*JV ' 2 ^V 



D (2) 



M 2 



Pn 



p-q 



Vv 



p-q 



(5) 



Performing the Wick rotation v — > iv and the variable 
transformation Q 2 = q 2 + i/ 2 , the self-energy becomes 



SAP =— dQ 



+ l,V^^ T H 



dv 



Q 2 



M 



5M ct {K) (6) 



where <5M C *(A) derives from counterterms required for 
renormalization [34| and the Lorentz contracted Comp- 
ton tensor is 



(7a) 



= -3Q 2 h(iv,Q z )+ 1 + 2— \ QH 2 {iv,Q 2 ) . (7b) 



T£ = ST^iu, Q*)+i\-—\ T 2 (iu, Q") , 



Or 



The scalar functions {Ti, t{) can be evaluated using a dis- 
persive analysis. It is known the (T±,ti) functions re- 
quire a subtracted dispersive analysis while the (T 2 ,i2) 
functions can be evaluated with an unsubtracted disper- 
sion relation [32| ■ In Ref. |28| , it was claimed the elastic 
contributions to t\ could be evaluated with an unsub- 
tracted dispersive analysis. However, performing an un- 
subtracted dispersive analysis of the elastic contributions 
to Eqs. ||7|) by inserting a complete set of elastic states 
into Eq. (|3J), leads to inconsistent results: 
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(8a) 

(8b) 



o 
with r e i = -Tjj2 ■ If both parameterizations of the clastic 

contribution were to satisfy unsubtracted dispersion rela- 
tions, the following positive-definite integral would have 
to vanish 



3a [°° G 2 E (Q 2 )+r el G 2 M (Q 2 ) 



2n 



l + T e l 



(9) 



Equating Eqs. (|4aj) and (|4b[) allows one to solve for Ti in 
terms of U and vice versa, and to demonstrate that if the 
clastic contributions to Ti(ti) do not satisfy an unsub- 
tracted dispersive analysis, then neither will the elastic 
contributions to ii(Ti). Eq. ([8b")) was used in Ref. HH 
and is often quoted as the elastic contribution to the nu- 
cleon self-energy. 



Starting from either Eqs. (|7a[) or ()7b[) . performing a 
subtracted dispersive analysis of (Ti,ti), and an unsub- 
tracted analysis of {T^ti), using a mass- independent 
renormalization scheme (dimensional regularization) one 
arrives at 133 



5M 1 = SM el + SM l 
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(10) 



with 
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(12) 



where t = v 2 /Q 2 , Fi(is,Q 2 ) are the standard nucleon 
structure functions and v t h = w^ + (m 2 + Q 2 )/2M; 



5MS „6 
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dQ 2 T!(0,g 2 ), (13) 



and 



<5M C 
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(14) 



where C\.i are Wilson coefficients determined from the 
operator product expansion of the counterterms 34 1. 



The UV divergence has been entirely cancelled by the 
counterterm and 5M ct is a remaining finite contribution 
with residual scale dependence. The scales Ao and Ai 
can be chosen arbitrarily provided their values are in the 
asymptotic scaling region. Restricting our attention to 
the isospin breaking contribution, with 25 = m d — m u 



<5M p c l n = 3aln 



e 2 u m u 



e 2 d m d 



SwMS 



( P \5(i 



■dd)\p) 

(15) 



with e u = 2/3 and e d = —1/3. In QCD, m Utd ~ 5, so the 
entire contribution is numerically second order in isospin 
breakin g, Q (a5), and for practical purposes can be ne- 
glected [34[ • Estimating the size of this term, with A 2 = 
100 GeV 2 , A 2 , = 2 GeV 2 yields \5Mf_J < 0.02 MeV. 

The remaining contribution to the self-energy is the 
subtraction term, which can not be directly related to 
experimentally measured cross sections. We now have a 
better theoretical understanding of this term enabling 
a more robust determination of its contribution than 
has been previously made. While the function is not 
known, the low and high Q 2 limits can be determined 
in a model independent fashion; the asymptotic region is 
constrained by the operator product expansion (OPE) to 
scale as liniQ_ i . 00 Ti(0, Q 2 ) ~ l/Q 2 1341 w hile the low Q 2 
limit is fixed by non-relativistic QED [35H39J 



t 1 (o,q 2 ) = 2k(2 + k)-q 2 !^[(i 



K ) r M - r l 



W 2 ~ 2M '^ 



(16) 



where k = -F2(0) is the anomalous magnetic moment, 
rE(rM) is proportional to the slope of the electric (mag- 
netic) form factor and commonly denoted as the nucleon 
electric (magnetic) charge radius and /3j\/ is the magnetic 
polarizability. 

A direct evaluation of Eq. (Tf"3"|) with Eq. (TIB]) diverges 
quadratically resulting in an uncontrolled uncertainty. 
However, the displayed Q 2 dependence is not of the form 
required by the OPE, so we are necessarily led to intro- 
duce model dependence. The first few terms in Eq. (Tl7j]) 
are recognized as the low-Q 2 expansion of elastic form 
factors and the magnetic polarizability term is the lead- 
ing inelastic contribution. In evaluating the elastic con- 
tributions to T^ v only the elastic u-spinors need be used 
in the dispersion relation. If one uses the full Feynman 
propagator in the full amplitude, a procedure known to 
be correct in the point- limit (as for the electron), and 
vertex functions with ordinary F\ and Fi form factor 
contributions, then the specific elastic terms of Eq. p^|) 
would arise [39J, |40(. This suggests a re-summation in 
which one uses the appropriate elastic form factors. The 
inelastic contribution can be multiplied by a dipole form 
factor (1 + Q 2 /mg)~ 2 , such that it has the correct asymp- 
totic limits as Q 2 — > 0, oo. The parameter m 2 , should be a 



typical hadronic scale and we will take m\ = 0.71 GeV 2 . 
The subtraction term is then approximated by two pieces 
which have the correct low and high Q 2 limiting behav- 



ior, 



Ti(0,Q 2 ) 



2G 2 M (Q 2 ) 

+ Q 2 2M^ 
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leading to the convenient separation 
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(18a) 
(18b) 



The second term, generated using the model assumptions 
described above, will cause the largest uncertainties, as 
we show below. 

Evaluation of contributions- In all subsequent evalua- 
tions, we take A 2 , = 2 GeV 2 for our central values and the 
range 1.5 2 < A 2 , < 2.5 GeV 2 to estimate uncertainties. 
We begin with an evaluation of the clastic contribution, 
Eq. (jlip . The form factors are well measured over the 
kinematic range required by the integrals, which are rep- 
resented by a number of analytic fits. The elastic contri- 
butions converge well at the upper limit, which may be 
taken to infinity with negligible error. Using the Kelly 
parameterization of the form factors [4l| , or an updated 
version [42T - |44| . the elastic contribution is given by 



5M e 



1.39(02) MeV. 



(19) 



The uncertainty is determined through an uncorrelatcd 
Monte-Carlo evaluation of the fit parameters in the 
parametrization.lt is also interesting to note, that if the 
simple dipole parameterization of the form factors is 
used, the same value within the quoted uncertainty is 
obtained. 

In the inelastic contribution, Eq. (TT2"|) . most of the 
support for the integrals lies in the resonance region, 
where there are good data from JLab, and there are 
analytic fits valid in the resonance region for both the 
neutron and proton structure functions from Bosted and 
Christy [45|, |46| (we also remind the reader the neu- 
tron functions are determined from deuterium-Compton 
scattering with the additional uncertainties captured in 
the coefficients of the neutron functions, and propagated 
into our uncertainties through a Monte-Carlo treatment). 
Their quoted range of validity includes Q 2 up to 8 GeV 2 
and W up to 3.1 GeV (W 2 = M 2 + 2Mv - Q 2 ). To 
extend the W range, we use the parameterizations of 
Refs. [47|, [48( which fit proton structure functions in the 
diffraction region using forms recognizable as Pomeron 
and rho meson Regge trajectories. The former is isoscalar 
and the latter isovector, so we have a straightforward ex- 
tension to the neutron case. Taking A 2 , = 2 GeV 2 and 



Wtrans = 3.1 GeV as the transition between the two pa- 
rameterizations, the inelastic contribution is 



SM 



inel I 



0.057(16) MeV. 



(20) 



The uncertainties are estimated by the range of Ag given 
above as well as varying the transition value of W be- 
tween 2.5 < Wtrans < 3.5 GeV. These two variations 
dominate the uncertainty estimate. The numerical inte- 
gration is insensitive to the upper limit of the W integra- 
tion through W max ~ 200 GeV (or x ~ 10~ 4 ). 

We are left with the subtraction terms. Using the 
model assumptions described above, the contribution 
from the elastic subtraction term, Eq. (|18a[) . is 



SM, 



sub I 



\p— n 



-0.62(02) MeV 



(21) 



It is interesting to note the sum of Eqs. (fT9)) and (f2"Tj) 
is surprisingly close to that of Ref. [28[ (although the 
individual proton and neutron elastic self energies are 
different). 

The most troublesome contribution to evaluate is that 
of the inelastic subtraction term, Eq. (|18b[) . This contri- 
bution is proportional to the isovector nucleon magnetic 
polarizability /3 p _„. The determination of this isovec- 
tor quantity was part of the motivation for the recent 
deuterium Compton scattering experiment, MAX-Lab at 
Lund |49f . for which we are still awaiting results. The 
HIGS experiment [50| at TUNL will also help determine 
this quantity. From chiral perturbation theory, one ex- 
pects the isovector polarizabilities to be small; the lead- 
ing contribution to the polarizabilities occurs at order 
P 3 and these are purely isoscalar. The isovector contri- 
butions arise at order P 4 and are suppressed in the chiral 
power counting [5l|. A recent review provides the con- 



servative estimate /?„ 



-l±lx 10~ 4 faV [30J. Using 



this in Eq. (|18b[) provides the determination 



SM? 



= 0.47 ± 0.47 MeV , 



(22) 



(a smaller value of ttjq would reduce these values) . 

Adding all the various contributions, Eqs. (|T§)) , ([2"0"|) . 
(|2Tj) and ([2"2"|). we arrive at 



SM 1 



p—n 



1.30(03)(47) MeV, 



(23) 



where the second uncertainty arises from the inelastic 
contribution to the subtraction term. Clearly, any im- 
provement in our knowledge of /3 p _ n will significantly 
improve our ability to determine the electromagnetic con- 
tribution to M p — M n . 

The isovector magnetic polarizability- Within the 
model assumptions used to arrive at Eqs. (TT8"|) . we can 
combine the experimental value for M n — M p with lat- 
tice QCD determinations of the m^ — m u contribu- 
tion. There are three published numbers from lattice 
QCD [1J, [l3, LL7J]: which are uncorrelated. For each re- 
sult, we combine the quoted uncertainties in quadrature, 
and then perform a simple weighted mean, arriving at 



SM. 



latt 

m d -m u |p_ rl 



-2.53(40) MeV. 



(24) 



Combining this with Eqs. <JTJ> , (fl8b)) . (fT9l) . (|20l) . (|2T1) and 
our value for m,Q, we find 



/3 p _„ = -0.87(85) x 10~ 4 fm 3 



(25) 



in good agreement with current estimates [30|. 

Model independence- One can infer the nucleon isovec- 
tor electromagnetic self-energy without recourse to mod- 
els by utilizing the known mass splitting, Eq. (fTJ, com- 
bined with the lattice QCD determination of the of the 
contribution from md — m u , Eq. (|24j) . 



(26) 



8Ml_ n = 1.24(40) MeV. 



Combined with Eqs. ([TTJl) and (|2U)) . this can be trans- 
lated into a model-independent bound on the unknown 
subtraction function 



3a 
IQ-kM 



dQ 2 Tf~"(0,Q 2 ) = 0.21(02)(40) MeV. (27) 



This is compared with Eqs. (j2"Tj) and (j2"2"j) which give 
0.15(02)(47) MeV for the same quantity. This bound 
demonstrates that our treatment of the subtraction func- 
tion, while not model-independent, is also not wildly 
speculative, but in agreement with the combined con- 
straint of experiment and lattice QCD. 

Conclusions- We have provided a modern and ro- 
bust determination of the isovector electromagnetic self- 
energy contribution, 8M^_ n = 1.30(03)(47). A technical 
oversight in the evaluation of the elastic contribution was 
highlighted resulting in a larger central value than pre- 
viously obtained [28|]. Modern knowledge of the struc- 
ture functions was used to constrain the elastic and in- 
elastic contributions, reducing the uncertainty from these 
sources by an order of magnitude (±0.30 MeV [28[ com- 
pared to our ±0.03 MeV). However, a careful analysis 
of the subtraction function has yielded an overall larger 
uncertainty than previously recognized. The larger cen- 
tral value suggests a larger contribution to M p — M n 
from m,d — m u , consistent with expectations from lattice 
QCD, thus impacting the phenomenology of Rcfs. |22| - 
|27J . With plausible model assumptions and additional 
input from lattice QCD, this knowledge can be used to 
provide a competitive estimate of the nucleon isovector 
magnetic polarizability, albeit still with a 100% uncer- 
tainty. Alternatively, a bound can be placed on the un- 
known subtraction function, which can not otherwise be 
determined, and lends further support for our determi- 
nation of /3 p _ n . 
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